Extracting the charm cross-section from 
semileptonic decays into muons 



Haidong Liu a,b ' c Yifei Zhang b ' c Chen Zhong a ' d Zhangbu Xu a 

& Brookhaven National Laboratory, Upton, New York 11973 
b University of Science & Technology of China, Anhui 230027, China 
c Lawrence Berkeley National Laboratory, Berkeley, California 94-720 
^Shanghai Institute of Applied Physics, Shanghai 201800, P.R. China 



Abstract 

We propose a sensitive measurement of the charm total cross-section at RHIC 
through muon identification from charm semileptonic decay at low transverse mo- 
mentum (pt)- This can test the binary-collision scaling (A^ n ) properties of the 
charm total cross-section and be used to study whether heavy- flavor quarks, which 
are used as a probe, are produced exclusively at the initial impact in hadron-hadron 
collisions. The effect of the charm semileptonic decay form factor on extracting the 
total charm cross section and on the shape of the lepton spectra are also discussed in 
detail. We conclude that lepton spectra from charmed hadron decays at transverse 
momentum px — 1.0 GeV/c are sensitive to the charmed hadron spectrum shape. 
Therefore, the interactions of heavy quarks with the medium created in relativistic 
heavy-ion collisions, especially the flow effects, can be extracted from the lepton 
spectra from charmed hadron decays at low px- 
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In relativistic heavy-ion collisions, charm quarks are believed to be produced 
in the early stage via initial gluon fusion and their production cross-section can 
be evaluated using perturbative QCD [1]. Study of the Nu n scaling properties 
of the charm total cross-section in p+p, d+Au and Au+Au collisions can test 
if heavy-flavor quarks, which are used as a probe, are produced exclusively at 
the initial impact. The interactions of heavy quarks with the medium provide 
a unique tool for probing the hot and dense matter created in ultra-relativistic 
heavy-ion collisions at the early times. At RHIC energies, heavy quark energy 
loss [2], charm quark coalescence [3,4,5,6], the effect of J/ij) production from 
charm quark coalescence on the interpretation of possible J/ip suppression 
due to color screening [7], and charm flow [8,9,10] have been proposed as 



Preprint submitted to Elsevier Science 



8 February 2008 



important tools in studying the properties of matter created in heavy ion 
collisions. The last three effects depend strongly on the charm total cross- 
section and spectrum at low px- 

However, it is difficult to directly reconstruct charmed hadrons and single elec- 
trons from charm semileptonic decay in hadron-hadron collisions with high 
precision at low pr, where the yield accounts for a large fraction of the total 
cross section [11,12,13]. The difficulties are due to large combinatorial back- 
grounds in charmed hadron decay channels, and the overwhelming photon 
conversions in the detector material, and n° Dalitz decays in electron de- 
tection. Nevertheless, the charm total cross-sections have been measured in 
d+Au collisions at RHIC by a combination of the directly reconstructed low 
p T D° — > Krr and the non-photonic electron spectra [11], and by electron 
spectra alone at p T > 0.8 GeV/c [13,14]. Although the systematic and sta- 
tistical errors are large, the result indicates a much larger charm yield than 
predicted by pQCD calculations [11,1]. Recently, the measurements of high p-p 
electrons from heavy-flavor semileptonic decays have posed challenges to our 
understanding of partonic energy loss in the medium [13,16,15,17]. Since most 
of the measurements to date at RHIC are from indirect heavy-flavor semilep- 
tonic decays, it is therefore important to find novel approaches to improve the 
measurements and also study in detail how to extract the maximum informa- 
tion about the heavy-flavor spectrum from its lepton spectrum. This includes 
studying the effects of the form factors of charmed hadron semileptonic decays 
on the lepton spectra [18]. 

In this paper, we propose a new method to extract the charm total cross- 
section by measuring muons from charmed hadron semileptonic decay at low 
Pt (e.g. 0.16^pr^0.26 GeV/c). Since muons in this pt range are a very uniform 
sample of the whole charmed hadron spectrum, the inferred charm total cross- 
section is insensitive to the detail of the charm spectrum. Once the cross- 
section is determined, the electron spectrum at higher p T can be used to 
sensitively infer the charmed hadron spectral shape. Meanwhile, we survey 
the form factors used in charm semileptonic decays generated from Particle 
Data Group [18], in the PYTHIA event generator [11,13,14,19], by pQCD 
predictions [1] and from the CLEO inclusive measurement [20]. We find that 
the lepton spectra from these different form factors can be different by a factor 
of 1.5. 

We first study the charmed hadron semileptonic decay form factor and its ef- 
fect on the lepton spectrum. Fig. 1 shows the electron momentum spectra from 
charmed meson decays at rest generated using the Particle Data Table [18], 
PYTHIA [11,13,14,19], pQCD calculations [1] and from the CLEO preliminary 
inclusive measurement [20]. The spectrum generated by the PDG is according 
to the form factor of charmed meson decays to pseudoscalar K + 1 + v, vector 
meson K* + / + v and non-resonance (Kir) + 1 + v where the K* mass is used for 
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the (Kir) system. The decay partial widths (T) of the three dominant decay 
channels are: 



(1) K + / + v with pseudo scalar meson in final state (D ± 7.8% 

dT p\ 

oc 



dq 2 (l-q 2 /M* 2 ) 2 
where q 2 is the invariant mass of the virtual W — > ul, p^ is the momentum 
of the kaon, and M* = 0.189 GeV/c 2 is a parameterization of the effective 
pole mass in the decay. 

(2) (Kir) + with non-resonant Kir in final state(D ± 4.0%) 

We use the K* mass (0.892 GeV/c 2 ) for the Kir invariant mass and the 
form factor is the same as in the decay to the pseudo scalar meson. 

(3) K* + I + v with vector meson in final state(D ± 5.5%) 
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where ^/ is the decay angle between the lepton and the vector meson, 
Pv is the vector meson momentum, 



H ± (q 2 ) = (M + m)A 1 (q 2 ) + ^f-V(q 2 ) 



and 



1 dM 2 n 2 

Ho(q 2 ) = 7T~l(M 2 - m 2 - g 2 )(M + m)A l (q 2 ) ~ i^Mf)] 
Iraq M + m 

where A 1<2 , V take the form of l/(l-q 2 /Mf v ) with M\ = 2.5 GeV/c 2 , 
My = 2.1 GeV/c 2 andr y = V(0)/A 1 (0) = 1.62±0.08,r 2 = A 2 (0)/A 1 (0) = 
0.83 ±0.05 [18,21]. 

For different charmed hadrons, we assume that the relative branching ratios 
among these three channels are the same, and their decay electron spectra 
are the same. The overall charmed hadron to electron branching ratio T(c — > 
e)/r(c — > anything) is 10.3% [18]. There is a possible ~ 5% difference between 
electron and muon decays due to phase space which was not taken into account 
in this analysis. Electrons at high momentum are mainly from decay channel 
(1) K + l + v because the kaon is lighter than the K* and the form factor of the 
decay channel to a vector meson (K* + 1 + v) favors a low momentum lepton 
and higher momentum neutrino. Since PYTHIA uses a simplified vector meson 
decay form factor [19], it tends to produce a softer electron spectrum. Both 
the parameterization by Cacciari [1] and formulae from the PDG agree with 
CLEO's preliminary electron spectrum. In addition, we also find that although 
the charmed mesons (D ± and D°) from ^(3770) decay have a momentum of 
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Fig. 1. Electron momentum spectra from charmed meson decays at rest. Histogram 
is the one with form factor from PDG. The dashed line is that of Cacciari's pa- 
rameterization. The dotted line is from the simplified vector-meson form factor in 
PYTHIA. The solid line is that from the PDG which takes a form factor from 
(3770)-> D -> e. 



244 MeV/c only and without correction of final state radiation [20], it affects 
slightly its subsequent electron spectrum. 

We use our electron momentum spectra and that of Cacciari to generate 
electron spectra from charmed decay at RHIC. A power-law function of the 
charmed hadron transverse momentum spectrum was obtained from minimum- 
bias Au+Au collisions [16] . Fig. 2 shows the ratios of those electron spectra di- 
vided by the spectrum using PYTHIA decay form factors [19,11]. The slightly 
soft form factor of the charm semileptonic decay in PYTHIA convoluted with 
a steeply falling charm spectrum produces an electron px spectrum in Au+Au 
collisions at RHIC, which can be lower than the correct one by up to a factor 
of 1.5 at high p?. Part of the discrepancy between experimental results and 
PYTHIA in electron spectra [11,13] can be explained by the decay form factor. 
Taking D ± and D° from ^(3770) decay as if it were at rest, results in slight 
change on the electron spectrum. 

Since the parameters for the form factor from the PDG are extracted from 
experimental data, and the form factor describes data well, we will select this 
form factor with D at rest, to generate lepton spectra for the remainder of 
the paper unless otherwise specified. A power-law function is used to create 
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Fig. 2. Charm-decay electron spectra for three different form factors divided by the 
spectrum using the PYTHIA form factor. The histogram represents the ratio using 
the PDG form factor. The solid line is Cacciari's parameterization, and the dashed 
line is from PDG and takes a form factor from ^(3770)— > D — > e. See text for detail. 

charmed hadron p? spectra. The function takes the form: 

dN _dN2{n-l){n-2) 2p T 
2ndyp T dpT dy n (n — 3) 2 (pt) 2 (PT){ n — 3) 

where dN/dy is the yield and n and (pt) are the parameters controlling the 
shape of the spectrum. Fig. 3 shows the charmed hadron p? spectrum be- 
fore and after requiring its decayed muons at 0.16 < pt < 0.26 GeV/c. The 
similarity of the spectral shape shows that the muon selection reasonably uni- 
formly samples the entire charmed hadron spectrum. The muons in this px 
range sample 14% of the charmed hadron spectrum. The muon yield is about 
1/70 of the charmed yield due to an additional 10.5% decay branching ratio. 
Fig. 4 shows the dependence of the muon yield on (px) for a fixed total charm 
yield. The yield is normalized to yields at n — 10 and (p?) — 1-3 GeV/c. We 
also note that the muon yield has a very weak dependence on n which demon- 
strates that over a wide range in (pr), the muon yield is within ±15%. This 
is in contrast to the large variation of the electron yield integrated above pt 
of 1.0 GeV/c, where a factor of 8 variation is seen in Fig. 4. This may explain 
the difference between results from Au+Au, p+p and d+Au collisions [22], 
pointing to a softer charm spectrum in Au+Au collisions. On the other hand, 
once dN/dy is determined by the low pt muons, the electron yield at higher 
Pt is very sensitive to the (px) of the charmed hadron spectrum. If the lepton 
spectrum at low px can be measured, the electron yield at higher p T essentially 
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Fig. 3. Charmed hadron (D) spectra {dN / dydpx) as a function of pt at midrapidity 
before (solid line) and after (dashed line) a muon selection of 0.16 < pt < 0.26 
GeV/c and \yi\ < 0.5. The later was scaled up by a factor of 70. 

determines the charm (pr)- In fact, the dependence on (p?) is almost linear. 
At lower beam energies at RHIC, the electron background is too high to make 
a meaningful charm measurement [23] due to the overwhelming photon con- 
versions in the detector material and background from 7r° Dalitz decays. The 
alternative muon measurement can be used to extract the charm total cross- 
section and study the energy excitation function of charm production [11,1]. 
In central Au+Au collisions, c — > fi in this pt range is estimated to be about 
0.1 per event due to the high production of charm quarks at RHIC. If the 
charm-decay muon yields are verified to follow binary collision scaling in p+p, 
d+Au, minimum-bias and central Au+Au collisions, we will be able to use 
it as an experimental measure of for peripheral collisions where Glauber 
model calculations [24] result in a large uncertainty on Nf, in . 

The key issue then is how to identify and measure the charm lepton spec- 
trum at very low px in the absence of high precision measurements of directly 
reconstructed charmed hadrons. First, we will demonstrate the feasibility by 
simple estimates and then using more detailed Monte Carlo simulations from 
a combination of two types of detectors: the Time Projection Chamber (TPC) 
and Time-of-Flight (TOF) detectors. 

The ionization energy loss difference between n and \i is about 12% at p = 0.2 
GeV/c and is much smaller at higher momentum [18]. Typical TPC dE/dx 
resolutions are 4% at NA49/SPS, 8% at STAR/RHIC and a projected 6% 
at ALICE/LHC. This means that the separation between tt and \i is about 
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Fig. 4. Lepton yields relative to the fixed total charm cross section as function of 
power-law parameters {px) for a charmed hadron transverse momentum spectrum. 
Solid line shows muon yields with a kinematics selection 0.16 < px < 0.26 GeVc 
and \yi\ < 0.5. Dashed line shows electron yields with px > 1.0 GeV/c. 

1.5cr to 3(7. Typical TOF timing resolution is about a t = lOOps at a distance 
of / = 2m. Therefore, the separation of tt and /j, measured in time is: At = 
{\jm\jp 2 + 1 — ^m 2 /p 2 + 1) x / x c ~ 600 ps. This is about a 6<r separation, 
which decreases rapidly as 7 3 at higher momentum. 

For our more detailed simulations, we used the HIJING event generator [25] 
to simulate Au+Au collisions in the STAR detector at RHIC to demonstrate 
how muons can be identified and how the background can be subtracted. The 
same approach can be used for detectors with similar configurations, such as 
CDF and ALICE. Muons were identified by measuring the energy loss in the 
Time Projection Chamber and the velocity in the Time-of-Flight patch at 
STAR [23,26]. From the dE/dx and Time-of-Flight resolutions, we conclude 
that the dE/dx difference between pions and muons is about 1.5(7 and the 
mass resolution is about 5<r at pt — 0.2 GeV/c. This provides a pion rejection 
factor of > 500, ~ 25 from each detector. The left panel of Fig. 5 shows 
the m 2 = (p/P/'y) 2 distribution from TOF after TPC dE/dx selections [23] 
for muon candidates and pure pion candidates. A clean muon peak can be 
identified within a mass window of 0.008 < m 2 < 0.014. The tail of the pion 
background can be evaluated by selecting pure pion candidates with a dE/dx 
cut and the result is shown as a dashed histogram in Fig. 5. The residual 
pions can be subtracted statistically from the distribution of the distance of 
the closest approach (DCA) to the collision vertex with the muon mass window 
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applying to the pure pion sample [26]. 



The STAR Collaboration have performed analyses of all the VO and other 
weak-decay strange hadrons and has published many papers on these top- 
ics [27]. Those decay topologies have shorter decay distances and comparisons 
between data and simulation show good agreement. Our analyses of DC A 
distributions are the similar and should lead to a good match between data 
and simulation [26]. Measurements of the muon spectra from charmed hadron 
semileptonic decays at low p T are not affected by the n° Dalitz decays and 
photon conversions that an equivalent electron measurement is. The dominant 
background muons from pion/kaon weak decays are subtracted using the DCA 
distribution. Other sources of background (p — > ,rj — > / ~ffi + fi~, etc.) are 
found to be negligible from the HIJING simulation. Fig. 5 shows the DCA 
distribution of muons from pion and kaon decays in minimum-bias Au+Au 
collisions. Also shown is the DCA distribution of particles at the same p T orig- 
inating from the collision vertex. The integrated yields of these two distribu- 
tions reflect the appropriate muon(from charm)-to-hadron ratio with hadron 
yields estimated from HIJING and the charm-to-muon yield from d+Au colli- 
sions scaled by the corresponding N bin . It is clear that the two distributions are 
very different and muons from charm decaycan be reliably obtained from this 
method with high precision when the DCA distribution from real data is fitted 
to the combination of these signal and background distributions. Since charm 
production is expected to scale with Nu n while low-py pions scale with number 
of the participant nucleons,iVp art , we expect a better signal-to-background ra- 
tio in central Au+Au collisions. Additional vertex detectors, such as the Silicon 
Strip Detector (SSD) [28] and/or the Heavy Flavor Tracker (HFT) [29], can 
greatly enhance this capability because the DCA of the muon from pion/kaon 
decay is determined by the decay distance and track geometry while for the 
muon from charm decay, the DCA accuracy is currently due to tracking resolu- 
tion and can be reduced to <100/im with a good vertex detector. In principle, 
the 20% systematic uncertainties presented by the STAR Collaboration [26] 
can probably be further reduced. It will then provide a complementary charm 
measurement in the era of directly reconstructing charmed hadrons from a 
displaced secondary vertex. 

With the improved charm semileptonic decay form factor [1] and possible 
high precision measurements of muons at low p?, we re-evaluate the sensitiv- 
ity of the lepton spectrum to the original charm spectrum as has been per- 
formed previously [10]. In our approach, we take the extracted charm spectrum 
from d+Au collisions as a baseline charm spectrum for nucleon-nucleon colli- 
sions [11]. A power-law function was fitted to the D° p T spectrum combined 
with its associated decay electron spectrum for the d+Au collision data [11]. 
A blast-wave similar to Ref. [10] with freeze-out temperature and freeze-out 
velocity parameters obtained from multistrange baryons {Tf = 160 MeV, 
Pmax = 0.6) was used to generate a charm spectrum [30]. The blast-wave as- 
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Fig. 5. Left panel: Particle mass m 2 = (p/f3/~/) 2 from Time-of-Flight detector af- 
ter TPC dE/dx selections of muon candidates (solid) and primary pion candidates 
(dashed). Right panel: Primary particle DC A distributon (dashed line) and muon 
DC A distribution from background after TPC dE/dx and TOF m? selections from 
HIJING simulations through realistic STAR detector configuration. The DCA dis- 
tritions of muons from charmed hadron decay are practically identical to the primary 
particle in this detector configuration. 



sumption is likely not appropriate for lepton > 2 GeV/c and therefore 
we focused our study for pt > 2 GeV/c. The spectrum was divided by the 
spectrum from d+Au. Fig. 6 shows the ratios of these spectra to the baseline 
spectra of charmed hadrons and leptons, respectively. The semileptonic decay 
greatly smears the spectrum and reduces the difference between the different 
spectrum shapes. However, it is clear that a reasonably realistic blast-wave pa- 
rameterization of charmed mesons in Au+Au collisions is very different from 
that in d+Au collisions. There is also a significant difference between spectra 
with different flow (blast-wave function) parameters at 0.5 < pt < 1.5 GeV/c. 
Fig. 6 shows that there is a factor of 3 difference at pt = 1.5 GeV/c between 
late freeze-out (T/ = 100 MeV, j3 max = 0.9) and early freeze-out (Tf a = 160 
MeV, P m ax = 0.6). Current measurements of non-photonic electron spectra 
and direct charm spectra seem to be consistent with a decreasing trend even 
at pr — 1.0 GeV/c [16], which is likely due to multiple collisions and ther- 
malization at low p T with early freeze-out [9] and not due to pQCD energy 
loss. However, since the overall normalization is not known, an early freeze- 
out scenario can be interpreted as suppression of the charm total cross-section 
as well. This ambiguity can be resolved by a measurement of the total cross- 
section. In addition, the errors on the current measurements are large in this 
Pt range. We advocate improving the measurements of electrons in this pt 
range to assess if charm thermalizes in the medium and has similar flow and 
freeze-out as multistrange hadrons [8]. 

In summary, we propose a sensitive measurement of the charm total cross- 
section at RHIC. The charm-decay muons at pt — 0.2 GeV/c can be identified 
by energy loss in a Time Projection Chamber and velocity from a Time-of- 
Flight detector. The background from pion/kaon decays can be subtracted 
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Fig. 6. Charmed meson spectrum of blast-wave form with Tt Q = 160 MeV and 
Pmax = 0.6 (dashed line) divided by the spectrum with a power-law parameteri- 
zation of d+Au data from STAR [11], and corresponding electron spectrum ratio. 
These represent an early freeze-out scenario. The dash-dotted line is for an electron 
spectrum from a later freeze-out assumption with Tf = 100 MeV and (3 max = 0.9. 

using the distribution of the distance of the closest approach to the collision 
vertex. We present simulations using HIJING within the realistic environment 
of the STAR detector. Detailed studies show that charm-decay muon spectra 
at low pt are proportional to the total charm cross-section. There is very weak 
dependence of the muon yields at the level of ±15% on charmed hadron shape 
over a wide range of spectrum parameters. Detailed comparison of different 
form factors of the charm semileptonic decay shows up to a factor of 1.5 
difference in the resulting lepton spectra between PYTHIA and more realistic 
form factors. A power-law spectrum and a blast-wave-type spectrum possibly 
present in Au+Au is different by about 30% at pt — 1 GeV/c. We conclude 
that the lepton spectrum at low pt from charmed hadron decays is sensitive 
to the charm spectrum. Therefore, electron spectra at 0.5 < pr < 1-5 GeV/c 
can be used to study charm quark thermalization and radial flow. 
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